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WITH  EXTENDED  RESONANT  MEDIUM 


I.  RESEARCH  OBJECTIVES 


In  1954  Dicke1  pointed  out  the  possibility  of  cooperative 

spontaneous  emission  of  radiation  from  an  excited  system  of  N  two- 

o 

level  atoms  at  a  rate  proportional  to  N  rather  than  N  which  is 

the  case  for  non-cooperative  emission.  This  phenomenon  is  called 

superradiance  or  superfluorescence  (SF).  It  can  be  utilized  to 

generate  coherent  short  pulses  of  high  intensity  from  far  infrared 

to,  perhaps,  x-ray  region.  It  has  been  under  very  intensive 

2 

theoretical  studies  since  1967  and  also  very  intensive  experimental 

3 

studies  since  1973.  These  recent  experimental  observations  have 
raised  some  serious  questions  about  old  theory  and  call  for  new 
approaches. 

The  general  purpose  of  this  research  project  has  been  to  try 
to  answer  some  of  these  questions  and  to  gain  better  understanding, 
through  analytical  studies,  the  following  aspects  of  SF: 


A.  Initial  Evolution: 

The  initial  stage  is  the  most  intriguing  and  most  interesting 
part  of  SF  for  the  following  reasons: 

1.  It  starts  as  incoherent  emission  and  spontaneously  evoles 
into  coherent  emission,  a  good  example  of  self-organization  which 

4 

is  the  subject  of  a  new  discipline  "Synergetics"  advocated  by  Haken. 
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Some  authors  have  gone  so  far  as  to  declare  that  the  way  this  tran¬ 
sition  is  made  is  the  only  interesting  feature  of  SF. 

2.  Superfluorescence  is  initiated  by  zero-point  fluctuations 
which  can  only  be  understood  quantum  mechanically.  But  in  the  semi- 
classical  approach,  the  role  of  this  fluctuation  is  simulated  by  an 
effective  initial  tipping  angle  8o  of  the  so-called  Bloch  vector. 
This  angle  has  become  a  hot  issue  of  current  debate  because  widely 
different  values  have  been  predicted. 

3.  In  experimental  observations,  the  shot-to-shot  fluctuations 
are  very  large.  It  is  believed  that  this  tremendous  fluctuations 
arise  from  the  random  initiation  rather  than  from  long-term  stochas¬ 
tic  development. 


B.  Overall  Behavior: 

All  the  previous  calculations  resulted  in  single-pulse  SF.  But 
the  very  first  experiment  by  Skrabanowitz  et  al.  in  1973  gave  multi¬ 
pulse  oscillating  emission.  This  so-called  ringing  is  due  to  the 
effects  of  radiation  reaction  ignored  previously.  Inclusion  of  these 
effects  change  the  governing  master  equation  from  Markovian  to  non- 
Markovian.  Satisfactory  solution  of  such  non-Markovian  master  equa¬ 
tion  is  a  formidable  task. 
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ACCOMPLISHMENTS  OF  THE  RESEARCH 


A.  Ineoherence-to-Coherence  Transition 

The  quantum  number  that  plays  the  central  role  in  the  descrip¬ 
tion  of  superradiance  is  the  so-called  cooperation  number  r.  When 
the  emitting  atoms  either  are  confined  to  a  region  smaller  than 
the  wavelength  of  the  emitted  radiation  or  interact  with  only  one 
radiation  mode,  r  is  conserved.  However,  in  the  interaction  of 
multimode  radiation  with  extended  medium,  r  may  or  may  not  be  con¬ 
served  in  the  process  of  emitting  a  photon.  So  we  have  both  r- 
conserving  (rC)  and  r-nonconserving  (rNC)  processes  to  consider. 
Radiation  emitted  through  rC  process  is  coherent  while  that  through 
rNC  process  is  incoherent. 

Assuming  the  excited  atomic  system  to  be  in  a  Dicke  state 
|r,m>,  we  can  calculate  the  emission  rates  through  rC  and  rNC 
processes,  respectively,  as  follows: 


T  (r,m)  =  (r+m) (r-m+1 )yy  =  (N-2  -n)(n+l)yy. 


TrNC(r,m)  -  (r+'°Kr+m-.lJY(1.1,) 

,  (»-a»— n)(N— 2t— »— 1)T(1— u),  (2) 

where  m  is  1/2  of  the  difference  of  the  number  of  atoms  in  the 
excited  state  and  the  ground  state,  y  is  the  initial  spontaneous 
emission  rate  of  a  single  excited  atom,  y  is  the  geometric  para¬ 
meter,  Z  and  n  are  the  numbers  of  photons  emitted  through  rNC  and 
rC  processes,  respectively;  namely, 


2 


5 


For  multimode  radiation  interacting  with  extended  system,  we 
have  typically  y  <<  1.  Under  such  condition,  we  can  see  from  Eqs. 
(1)  and  (2)  that  at  the  beginning,  when  both  i l  and  n  are  very 
small,  photons  are  emitted  predominantly  through  rNC  processes  and 
undirectionally .  Later  op,  as  n  increases,  rNC  processes,  which 
emit  photons  directionally,  become  more  and  more  important.  There¬ 
for,  there  is  a  transition  from  incoherence  to  coherence  in  the 
spontaneous  emission  which  occurs  when  both  H  and  n  are  still  much 
smaller  than  N.  This  transition  is  the  first  subject  of  our  inves¬ 


tigation. 


1.  Numerical  Study 


Let  P(£,n)  be  the  probabilty  that  when  Jl+n  photons  have  been 
emitted,  £  of  them  are  through  rNC  processes  and  n  of  them  through 
rC  processes.  We  can  write  down  a  partial  difference  equation  for 
P(£,n)  as  follows: 

PC*.”)  -  iSjk  ♦  rrdbrTo  p**-1-”)  ‘4> 

where  a  =  y/(l-y).  We  have  obtain  the  solution  to  Eq.  (4): 

P( £ , n )  =jj  n  (l+ia)_6ij  £  II  (l  +  ia)"*jn!  an  *  (5) 

where  8.  =0,  1,  2,  ...,  and  the  summation  is  over  all  possible 
1  n+1 

choices  of  8.  such  that  J  8-  -  £. 

1  i=l 

Equation  (5)  has  been  evaluated  numerically  by  using  the  so- 
called  steepest-descent  method.  The  result  has  been  reported  in 
Publication  No.  1. 
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2.  Analytical  Study 

Let  P0  (t)  be  the  probability  that,  at  time  t.  Si  photons 

have  been  emitted  through  rNC  processes  and  n  photons  through  rC 

processes.  The  differential  equation  for  P.  (t)  is 

x>  j  n 

£Pl,n(t)  *  +  NnYgP,.n-1^t) 

-  (NY(l-p)  +  N(n+l)Yy}P£>n(t)  (6) 

We  have  obtained  the  solution  to  Eq.  (6)  as  follows: 


P£  n(t)  =  e~NYt{l-e"NYyt}n{NY(l-y)t}h'. 


(7) 


From  this  solution  we  have  calculated  the  expectation  values 
of  the  photon  numbers 

<l>  =  NY(l-y)t,  <n>  -  eNyyt  -  1  (8) 

and  the  emission  intensities  in  terms  of  photon  number 


IrNC  =  NY(l-y), 


rrC 


=  Nyye 


Nyy  t 


(9) 


Let  the  transition  point  be  determined  by  the  condition 
IrNC  =  IrC;  then  we  have 


t  =  1  *n{i^} 
c  Nyy  y 


(10) 


The  result  of  this  analytical  study  has  been  reported  in 


Publication  No.  2. 
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B.  Exact  Solution  of  Superradiance  Master  Equation 

The  overall  behavior  is  predominantly  determined  by  rC  pro¬ 
cesses.  If  the  active  region  is  small  enough  that  we  can  assume 
that  a  photon  escape  out  of  the  system  as  soon  as  it  is  emitted. 
Under  these  conditions,  the  evolution  is  described  by  the  so-called 
superradiance  master  equation  (SME): 


^P(n,t)  =  (N-n+1)  nP(n-l.t)  -  (N-n )  (n+1 )  P(n ,  t ) ,  (11) 

n=0,  1,  2,  ...,  N  (=2r ) 

where  P(n,t)  is  the  probability  that  the  atomic  system  is  in  the 
Dicke  state  |r,m>  such  that  r-m  =  n.  We  have  obtained  the  exact 
solutions  to  Eq.  (11)  under  various  initial  conditions. 


1.  Complete  Initial  Excitation 


In  this  case  we  have  the  initial  condition  P(n,0)  =  and 

the  exact  solution  to  the  SME  has  been  found  to  be: 

n 


P(n<N/2 , t )  =  [  A  e-(N-i)(i+l)t 

i=0  ’ 

P(N/2<n<N,t)  =  N  l  2A  e~(N-i)(i+l)t 
i=0  n' 


(12) 


+  N/r  (B  .  +  C  .t}e-(N-i)(i+l)t 

i=N-n-l  n’  n>1 


N/ 2-1 


P(N,t)  =  1  +  JQ  (BN>i  +  CN>it}e“<N-i)<i+1>t 


where  we  have  defined 
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B 


n,i  = 


(-pN-H-1  Nt  n  i 


(N-n) ! (n-i) !i ! (N-i-1 ) ! ( i+n-N+1 ) ! 


r  N-i-l 

:  2  -  (N-21-1)  l 


k=i+l 


k( 


_1 _ .1 

n-k+1 ) 


n,i  = 


(-1 )N~n~1(N-2i-l )2  N!  n! 


(14) 


(15) 


(N-n ) ! (n-i ) ! i ! (N-i-1 ) ! (i+n-N+1 ) ! 

From  this  solution  we  have  calculated  various  expectation 
values  of  physical  quantities.  The  most  important  one  is  that  of 
the  photon  number 


<n>  _  _  N/2~X  (N-l)l  -(N-i ) ( i+1  )t 

N  “  L  (N-i-1)  !i!  e 

r  1  2  ~ l 

x  2  -  (N-21-1)  l  i  +  (N-2i-l )  t  I  (16) 

L  k=i+l  K  J 

For  practical  purpose,  we  have  also  obtained  the  following  appro¬ 
ximate  expression 

=  eNt/N  -  2!(eNt/N)2  +  3!(eNt/N)3  +  . 


■r. 


xe 


-x 


x+Ne 


-Nt 


dx 


(17) 


As  by-products,  we  have  also  discovered  many  identities  invol¬ 
ving  binomial  coefficients  and  Stirling  numbers  of  the  second  kind. 
These  identities  are  new,  as  far  as  we  know,  and  might  be  of  pure 
mathematical  interest. 


These  results  have  been  reported  in  Publication  No.  3. 


2.  Arbitrary  Initial  Excitation 


The  evolution  of  the  probability  distribution  of  the  atomic 
system  is  considered  as  a  Markov  process  described  by  a  time- 
evolution  matrix  (TEM).  The  exact  analytic  expressions  for  the 
elements  of  this  TEM  have  been  obtained  by  solving  Eq.  (11).  It  has 
been  proved  that  the  TEM  satisfies  the  normalization  condition, 
the  initial  condition,  and  the  Chapman-Kolmogorov  equation  as  re¬ 
quired.  The  exact  general  formulas  for  the  expectation  values  of 
various  physical  quantities,  expressed  in  terms  of  an  arbitrry 
initial  distribution,  have  been  obtained.  For  example,  the  emission 
intensity  as  a  function  of  time  has  been  found  to  be 


I  = 


r-1 

I 

i=0 


^e-(N-i)(i+l)t  |  (_1)1(N_i)(i+1)(N_2i_i) 


x  (*-i-i)J. 

*=i+l  £  Jl!(N-*-i-l)! 


.  y  P„(0)(-l)£(N-£)! 

! (N-i-i-l ) ! (i-1 ) ! 


f  -  N-i-1  , 

I (N-i) ( i+1 ) (2-(N-2i-l )  J 


k=i+l 


+(N-2i-l)^{(N-i)(i+l)t 


•>]| 


(18) 


where  P^O)  is  the  initial  probability  distribution. 

The  general  formulas  have  also  been  applied  to  the  specific  case 
of  superradiant  initial  state.  The  most  noteworthy  result  is  the 
normally  ordered  intensity  fluctuation 


°n(I)  =  |tftanh(x/2)  -  2}  tanh(x/2) 

These  results  have  been  reported  in  Publication  No.  4. 

[ 


(19) 


C.  Oscillatory  Superfluorescence 

For  an  extended  system,  an  emitted  photon  cannot  escape  out 
of  the  system  immediately.  So  photon  reaction  must  be  considered 
and  the  radiation  process  becomes  non-Markovian .  The  basic  equa¬ 
tion  for  such  oscillatory  superfluorescence  and  superradiance  is: 

=  ApW(t)  -  ige"t/2T2  [(aR+  +  a+R~) ,  eV^KO)] 

-  aV),  eAF(t-t’)[(aR++a+E-),  W(f)]J 

x  e_(t+t' )/2T 2  dt  ’  (20) 

where  (a)  W(t)  is  the  density  matrix  of  the  system  composed  of  N 
two-level  atoms  and  a  single  axial  mode  of  the  internal  field  (i.e., 
the  field  inside  the  active  volume);  (b)  a+  and  a  are  the  photon 

■f 

creation  and  annihilation  operators;  (c)  R  and  R  are  the  collec¬ 
tive  atomic  excitation  and  deexcitation  operators;  (d)  A^W(t)  = 
k|  [aW(t ) ,  a^]  +  H.c.j  whre  K  =  c/2L  (L  being  the  length  of  the  active 
region)  is  the  damping  constant,  so  Ap  describes  the  irreversible 
escape  of  photons  from  the  active  region;  (e)  g  is  the  atom-field 
coupling  constant  and  is  the  smallest  atomic  relaxation  time 
which  will  be  assumed  to  be  infinite  in  the  time  scale  of  our  problem. 

All  order  of  atom-field  correlations  have  been  included  in  Eq. 
(20).  A  complete  solution  of  this  equation  seems  to  be  impossible. 

We  have  obtained  solutions  that  include  zeroth  order  and  first 
order  correlations.  From  these  solutions  we  can  calculate  the  expec- 


tion  values  of  various  physical  observables.  For  example,  the  radi¬ 
ation  intensity  as  a  function  of  time  have  been  found  to  be 
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Vt) 


N/ 2-1 

l 

■  i=0 


[N-i-1  ,-i 
2  -  (N-21-1 )  s]Xo(i't) 

[xo(i,t)  -  (N-i)(i+l)Yo(i,t)]  J 


(N 

-  (N-2i-l ) 

N! 

x  (N-i-1) ! i ! 

for  the  case  of  zeroth  order  correlation;  where 


(21) 


XQ(i,t)  =  e~Kt/' 


£cos  w^t  +  sin  (22) 

YQ(i,t)  =  (K/2u>i)2I  e“Kt/2  [(1  +  2toit/K)  sin  wit  -  w;.t  cos  u)it] 

(23) 

and  where 


k  =  £fCI (N~i ) ( i+1 )  -  K2/4 


(24) 


To  obtain  the  radiation  intensity  I^Ct)  for  the  case  of  first 
order  correlation,  we  need  only  to  replace  the  XQ(i,t)  and  YQ(i,t) 
in  Eq.  (21)  by  X^(i,t)  and  Y1(i,t)  which  are  the  inverse  Laplace 
transforms  of  X-^i.s)  and  Y^(i,s),  respectively: 

X1(i,s)  =  KI/ js(s+K)(s+2K)  +  KI(N-i )(i+l)(3s+2K)}  (25) 

Y1(i,s)  =  (3s+2K)(X1(i,s)}2  (26) 

The  details  of  these  results  have  been  reported  in  Publication 
No.  5. 
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D.  Superradiance  in  a  Cavity 

With  possible  applications  to  the  problem  of  optical  bista- 
4 

bility  in  mind,  we  study  the  cooperative  effects  in  the  interaction 
of  matter  with  radiation  in  a  cavity  which  is  partially  open,  re¬ 
sulting  in  the  broadening  and  frequency  shift  of  cavity  modes. 

1.  Frequency  Shift  and  Decay  Time  of  Radiation  in  a  Cavity 

We  consider  a  small  Fabry-Perot  cavity  imbedded  in  a  larger 

one  with  perfectly  reflecting  walls,  receding  to  infinity  eventually 

to  represent  the  universe.  This  is  a  generalization  of  the  Lang- 

5 

Scully-Lamb  model  for  single-quasimode  laser  operation  in  the  sense 
that  both  mirrors,  instead  of  just  one,  of  the  inside  cavity  are 
semitransparent.  We  have  obtained  a  complete  set  of  orthonormal 
modes  of  the  universe.  Using  these  mode  functions,  we  have  studied 
the  leakage  of  a  monochromatic  radiation  out  of  the  inside  cavity. 

We  have  found  that  this  leakage  leads  to  a  frequency  shift 


a  decay  time 


(27) 


t  =  -25,/c  In  R 


(28) 


and  an  enhancement  factor,  defined  as  the  ratio  of  the  amplitude 
of  the  mode  function  inside  the  cavity  to  that  of  outside, 


=  /( 1+3R) / ( 1-R ) 


(29) 


where  21  is  the  length  of  the  cavity  and  R  is  the  reflectance  of 


the  mirrors. 
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The  details  of  the  derivation  have  been  reported  in  Publication 
No .  6 . 

2.  Transition  Probabilities  between  Mixed  Dicke  States 

The  cooperative  phenomena  in  the  interaction  of  radiation  with 
matter  are  usually  discussed  in  terms  of  the  so-called  Dicke  state 
jr,m,a>  which  is  a  pure  state  in  the  terminology  of  density  matrix. 
Such  a  pure  state  is  specified  by  a  degeneracy  index  a  in  addition 
to  the  angular  momentum  eigenvalues  r  and  m.  Since  the  physical 
meaning  of  a  is  quite  obscure,  it  is  more  realistic  to  consider  a 
mixed  Dicke  state  |r,m>  to  which  each  |r,m,a>  contributes • with 
equal  weight  and  random  phase.  Using  perturbation  approach  and  the 
diagrammatic  technique  developed  by  us,  we  calculated  the  various 
transition  probabilities  in  the  interaction  of  multimode  radiation 
with  resonant  atoms  in  a  Fabry-Perot  cavity  of  dimensions  much 
larger  than  the  wavelength  of  the  radiation. 

a.  r-Conserving  Processes: 

P(|r,m>  -*■  |r,m-l>)  =  Ayf5^(r+m)(r-m+l)  (30) 

where  A  is  the  Einstein  coefficient  for  spontaneous  emission,  y  is 
the  geometric  parameter,  and  3  is  the  enhancement  factor  given  in 
Eq.  (29). 

b.  r-Increasing  Processes: 

P(|r,m>  -  | r+1 ,m-l> )  =  A( 1-y ) ( r-m+1 ) ( r-m+2 ) (N/2-r )/4r2  (31) 

c.  r-Decreasing  Processes: 

P(|r,m>  -*•  |r-l,m-l>)  =  A(  1-y ) (r+m) ( r+m-1 ) (N/2-r) /4r2  (32) 

The  preliminary  results  of  this  calculation  have  been  reported 
at  the  APS  meeting  in  Washington,  D.C.,  24-27  April  1980. 
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E.  Initiation  of  Superfluorescence  during  Pumping 

In  all  the  theoretical  sudies  on  superf luorescence ,  the  compl¬ 
ete  initial  excitation  is  assumed  to  be  achieved  instantaneously. 

But  this  is  impossible  realistically.  So  we  have  proposed  the 
following  master  equation  to  describe  the  initiation  of  superfluo¬ 
rescence  which  allows  finite  duration  for  pumping: 

^■P(N,n,t)  =  pP(N-l,n,t)  -pP(N,n,t)  +  yy (N-n+1  )nP(N, n-1 , t ) 
-  yy(N-n )(n+l)P(N,n,t)  (33) 


where  P(N,n,t)  is  the  probability  that  N  atoms  have  been  excited 
and  n  photons  have  been  emitted  at  time  t.  The  first  two  terms  on 
the  right  hand  side  of  Eq.  (33)  describe  the  pumping  as  a  Poisson 
process  at  a  constant  rate  p;  and  the  last  two  terms  describe  the 
superfluorescence  with  varying  cooperation  number,  where  y  is  the 
geometric  parameter  and  y  is  the  single  atom  emission  rate. 

Under  the  condition  that  n  <<  N,  we  have  found  the  solution  to 


Eq.  (33)  to  be 


P(N.n.t) 


(p/UY)N  -PYt 
N! 


n  i 

l  (-D1 

i=0 


e-(i+l)yyt 

T+I 


34) 


From  this  solution,  we  have  calculated  the  expectation  value  of 
the  initial  tipping  angle  to  be 


9q  =  (2/N)^exp(Tp/4tR) 


(35) 


with  quantum  fluctuation  in  the  range 


9o”  =  9 0 { 1  ±  (1  -  e  Tp/2TR)}i 


(36) 
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where  x^  is  the  pumping  duration  and  xD  =  1/(uyN)  is  the  super- 
P  k 

radiance  time. 

6 

Using  the  experimental  conditions  of  Vrehen  and  Schuurmans: 

N  =  2xl08,  Tp  =  2  ns,  x^  -  0.4  ns,  we  can  calculate  0Q  -  3.5  x  10-4 
±  -4  -4 

and  0Q  =  0.7x10  to  4.9x10  ;  which  compare  favorably  with 

6  —4  +  —4  —3 

the  measured  values:  0^  =  5x10  and  0  =10  to  2.5x10 

o  o 

The  details  of  the  derivation  have  been  reported  in  Publication 
No.  7. 
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